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Abstract 

High order perturbation theory corrections to the superfluid transition temper- 
ature in a weakly interacting Fermi gas with repulsive interaction are calculated. 
This involves calculating the contributions of third and fourth order diagrams in the 
gas parameter app and taking into account effects of retardation. The contributions 
from both second, third and fourth orders result in the effective attraction in p-wave 
channel. It is shown that the critical temperature is mainly determined by second 
and third orders terms of perturbation theory. The dependence of the critical tem- 
perature on an external magnetic field is found. We discuss possible applications of 
the results to the diluted '^He-'^He mixtures and trapped neutral-atom Fermi gases. 

PACS: 67.60.-g, 74.20.-z, 74.20.Mn 

1 Introduction 

Unconventional mechanisms of Cooper pairing have recently started to attract greater 
attention. This is primarily related to the discovery of high-temperature superconducting 
(HTSC) systems, superconductivity in organic materials and heavy-fermion compounds, 
and also because of the search for superfluidity in of ^He-^He mixtures and in trapped 
atomic Fermi gases. However, HTSC and heavy-fermion systems belong to a class of 
strongly correlated systems whose theoretical analysis requires the development of new 
methods. At the same time, ^He-^He mixtures and trapped atomic Fermi gases can be 
described using the model of a weakly interacting Fermi gas. In this case the interparticle 
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interaction can be either attractive or repulsive. In the attractive case conventional singlet 
Cooper pairing takes place where the orbital momentum of the pair is £ = 0, for which the 
critical temperature was first calculated by Gor'kov and Melik-Barkhudarov In systems 
with repulsive interaction, the formation of ^ = Cooper pairs is clearly impossible and in 
order to investigate the existence of superfluidity, we need to study the possibility of £ 7^ 
Cooper pairing. The possible existence of superfluidity in Fermi systems with repulsion 
was first indicated by Kohn and Luttinger in 1965. In they examined the contribution 
of collective effects to the scattering amplitude in a particle-particle channel which lead 
to effective quasiparticle interaction at the Fermi surface via polarization of the Fermionic 
background. A principal role in the formation of attractive harmonics in the effective 
interaction and consequently the superfluidity is played by the Kohn singularity in the 
effective interaction. In the three-dimensional case, it has the form 

+ rreg(g'). (1) 

In coordinate space the Kohn singularity leads to alternating sign and magnitude dependent 
RKKY-type interaction between quasiparticles: 

reff^(r) ~ ^ cos(2pi.r + Lp). 

It should be noted that the above contribution to the effective interaction decreases over 
large distances more slowly that the bare interaction f/o(r — r') and consequently corre- 
sponds to the main contribution to the scattering amplitude in the limit of large momenta 

A simple extrapolation made by the authors of yields extremely low estimates for 
superfluid transition temperatures in the limit I — > 2: 10~^^ and 10~^^ for ^He and the 
electron subsystem in the metal, respectively. It was subsequently shown in P, that 
effective attraction occurs also for the angular momentum ^ = 1 which gives the following 
expression for the critical triplet-pairing temperature in the second order of perturbation 
theory: 

rr. ~ { Svr^ 1 . r 13.01 

^ ^^"P r 4(21n2-l)(ap.p l " ^^^^ T^l " 

where A = {2apF)/'^ is the gas parameter, a is the s-wave scattering length, pp is the 
Fermi momentum, and e is the energy parameter, of the order of the Fermi energy, which 
provides a cutoff at high energies. Substituting experimental values for ^He where triplet 
pairing takes place, gives good agreement with experiment: Td ~ lO-^K. (Obviously, the 
bare interaction in real '^He is far more complex than that in the considered model). The 
aim of the present paper is to determine the critical superfluid transition temperature of a 
weakly non-ideal Fermi gas with repulsive interparticle interaction up to the preexponential 
factor. For this purpose we calculate the irreducible vertex in the Cooper channel in the 
third and fourth orders of perturbation theory with respect to the gas parameter A. We 
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also allow for renormalization of the singular parts of the Green's functions (corrections 
associated with the Z-factor and the effective mass) in the Bethe-Salpeter equation (|) 
and take into account retardation effects (the influence of the frequency and momentum 
dependencies of the irreducible vertex). This article is organized as follows. In Section 
2 we derive and analyze an equation for the critical temperature in a weakly interacting 
Fermi gas with repulsion. In Section 3 we calculate the irreducible vertex in the Cooper 
channel in the second, third, and fourth orders of perturbation theory. In Section 4 we 
examine the contribution of retardation effects. In Sections 5 and 6 we give the final 
formula for the critical temperature and discuss the contribution of the bare scattering in 
the p-wave channel. In Section 7 we note the possibility of a strong enhancement of Td in 
an external magnetic field. In Section 8 we discuss possible experimental applications of 
the obtained results. In particular, we discuss the possibihty of triplet Cooper pairing in 
^He-^He mixtures and in a trapped neutral-atom Fermi gases at ultralow temperatures. 



2 Superfluid transition in a Fermi gas with repulsion 

We consider a weakly interacting Fermi gas described by the Hamiltonian 

H = Hq + Hint = X! (^P ~ ^pa^pa + 
a p 

+ 2 ^ O-paCtp'/jOp'+q/jap-qa, 

"/3pp'q 

where the indices a, /3 = 1, 2 label the system components which we assume to have equal 
masses m and concentrations ni^2 = p^/Qn"^ , /i is the chemical potential, and the constant 
g characterizes the interparticle interaction which we shall assume to be point-like (here 
and subsequently we put ^ = 1). The specific physical meaning of the two components 
depends on the particular system. For example, for ^He-^He mixtures it corresponds to an 
"upward" and "downward" projection of the spin, whereas in the case of trapped atomic 
gas, it corresponds to a hyperfine-structure component (or projection of the nuclear spin). 
The considered form of the interparticle interaction assumes that only s-wave scattering 
takes place in the system, characterized by the scattering length a. (In the leading order 
of perturbation theory a = mg/A7c.) The corresponding small dimensionless parameter, 
the gas parameter A, is given by A = 2\a\pF/7c. We subsequently show how the final result 
is modified in the presence of scattering in channels with nonzero orbital momenta. It is 
well known that the appearance of superfluid pairing is associated with the presence of a 
pole in the two-particle vertex function F in the particle-particle channel (Cooper channel) 
for zero total momentum and frequency [^j. This vertex function F is a solution of the 
Bethe-Salpeter integral equation (Fig. |I]): 

r(j9i, -pi;P3, -P3) = f (j9i, -pi;P3, -P3) - (3) 



-T J2 /r(pi,-pi;g, -g)G'(o;„,q)G'(-a;„,-q)F(g, -g;p3,-p3) 



(27r)3' 
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Figure 1: Bethe-Salpeter equation for complete vertex T. 



where F is the irreducible vertex in the Cooper channel (having no singularities at zero total 
momentum and frequency), G is the single-particle Green's function, and the arguments of 
the vertex functions denote the corresponding sets of Matsubara frequencies and momenta: 
q = (u;„,q), pi = (co'„i,pi), and so on. Note that in formula (|^) (and in the following 
formulas) we do not explicitly indicate the indices distinguishing the components of the 
Fermi gas (for example, F should be considered as Tais-ys, and so on). Writing them in 
explicit form does not present any difficulties. We also note that the nonsymmetrized (in 
terms of the component indices) irreducible vertex function is used in equation (^). The 
vertex functions F and F appearing in (^) are in fact functions of the Matsubara frequencies, 
the absolute values of the incoming and outgoing momenta, and the angle between them. 
For example, we have 

r(pi, -Pi;P3, -Ps) = T{uJi,UJ3, |pi|, |p3|,cos(6'p^p3)). 

Thus, expanding G and as a series in terms of Legendre polynomials 

f (. . . , cos(^^)) = ^(2/ + l)f ,(. . .)Pi{cos9), (4) 

1=0 

F(. . . , cos(^^)) = ^(2/ + 1)FK. . .)Pi{cos9) 

1=0 

and integrating over angles, we easily obtain from (|^) the following equation for the singular 
part F;*'' of the £-th harmonic of the vertex function: 

rS'H^i,^^3, |pi|, |p3|) = (5) 
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i [ -i^A^'^li^U^n,\Pl\,\q.\)G{Un,q)G{-U^ 

n=—oo \ J 

As usual, the critical temperature corresponds to the appearance of a nontrivial solution 
of this equation which is related to singular (logarithmic) behavior of the Cooper loop near 
the Fermi surface. Thus, in order to determine the critical temperature in the leading order 
in A, it is sufficient to set the frequencies to zero and the absolute values of the momenta 
to pf in all the vertex functions contained in (j^). We then have 

- ^ E / (Ifs ^(^«' ^)G{-^n, -q) - "^ZW In ^, (6) 

where vp = mp^/27r^ is the density of states at the Fermi surface, m* is the effective mass, 
and Z is the residue in the singular part of the single-particle Green's function. In equation 
(^) £ ~ is the cutoff parameter which depends on the behavior of F at high momenta 
and frequencies. Now equation (^) can be rewritten in the form 



Fr = F,.Z^ — ^^In - -Fr, (7) 

m \T/ 

where f; = fi{pi = ppi^i = 0) F^'*'' = Tf\pi = ppjOOi = 0), so that a nontrivial solution 
is only possible for F < and occurs at temperature T = Td where 

{1 m 1 , , 

For the case of a Fermi gas with interparticle attraction we have Fq — 47ra/m < 0, and 
the system is unstable with respect to traditional s-wave pairing {/el = 0). The superfluid 
transition temperature in this case was obtained in to within terms 0(A°) and is given 
by 

1 cf2V^^ { TT 1 r 1 



T.o^-e-y .,exp|-^|.0.28..exp|--|, (9) 

where C = 0.58 ... is the Euler constant. This expression only differs from the correspond- 
ing expression for Tco in BCS theory in that the Debye frequency ujd is replaced by ^i? in 
the preexponential factor. This replacement means that in our case the entire Fermi sphere 
and not only its vicinity of the order of the Debye frequency, is involved in the pairing. 
Note that in order to find the preexponential factor in we need to calculate Fq to within 
terms of the second order of perturbation theory. 

For repulsive interaction, a > 0, equation (^) for i = only has a trivial solution and 
s-wave pairing is impossible. In this case, superfiuid pairing will take place in the channel 
having the orbital momentum i for which F/ is negative and has the maximum absolute 
value. As it is well known [Q, the scattering amplitude of slow particles with the orbital 
momentum i for the short-range potential has the order of magnitude a{apY\ where p is 
the particle momentum and a is the s-wave scattering length. Thus, in our particular case 



5 



Pi p3 Pi p3 




* 1 » > ^ ^ ► 



-Pl -P3 -Pl -P3 

Figure 2: Second-order diagrams in terms of gas parameter for irreducible interaction. 



the corresponding contribution to F; with i > from scattering on the bare interparticle 
potential has a maximum for i = 1 and has the order (app)^ ~ A^. At the same time, many- 
particle effects associated with polarization processes of the Fermi background in a Fermi 
gas have the order and, therefore, dominates for / > 1 [Q. Corresponding diagrams for 
F in the second order with respect to A are plotted in Fig. J^. For our particular case of 
point interaction the first three diagrams cancel out so that F is completely determined by 
the last exchange diagram and is given by (we assumed T = 0) 

f(u;i, UJ3, Pl, ps) = ^ n(u;i + us, pi + ps), 
Km/ 

where 

In this expression n{p) is the Fermi particle distribution function for T = 0, ^{p) = 

From formula ([T0|) we can easily obtain an expression for the irreducible vertex at zero 
external frequencies and momenta lying on the Fermi surface. In terms of the angle 6 
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between pi and ps we have 



i/irf(0,0,pi =Pf,P3 =Pf,cos9) = ^Ii-^ 



v^(l + cos^), Vl-cos^ 
1 H , In 



— cos^^ 

As a result of the integration with Legendre polynomials we obtain 

z/pfi = -A^fl - 21n2) < 0. (11) 
5 

All the other partial components Ti with i > 1 also correspond to attraction but are smaller 
than Fi and their absolute values decrease rapidly with increasing i (see Thus, we 
conclude that a weakly interacting Fermi gas with interparticle repulsion is unstable with 
respect to triplet p-wave pairing. The corresponding critical temperature in the leading 
order with respect to A is 

= '""P {"^} = ''"^ {" (21n2-l) ^} • ^^2) 

It can be seen from this formula that in order to determine the preexponential factor 
e in equation (|^) we need to retain terms up to order A^. (This follows from the fact that 
since Fi begins from terms A^, to obtain terms of the order A*^ in the exponent, we need 
to know Fi up to terms A'^ includingly.) 

Note that the contribution of triple collisions can be neglected within the considered 
accuracy since it has the order A^ §6]. 



3 Contribution of higher orders of perturbation the- 
ory 

The irreducible vertex F in the third and fourth orders of perturbation theory is given 
by the diagrams shown in Figs. ^ and ^ respectively. The points on these diagrams 
correspond to antisymmetrized two-particle interaction. In expanded notation when the 
interaction is represented as a dashed line (as in Fig. ^) these corresponds to two different 
ways of connecting the incoming and outgoing lines. 

Figures ^ and ^ only give skeleton diagrams (without the self-energy insertions) and Fig. 
^ only gives " nonoriented" diagrams. The corresponding Feynman diagrams are obtained 
by arranging the arrows on the lines (taking into account the particle number conservation 
at the vertexes) and also the incoming and outgoing momenta. Figure ^ shows an example 
of such an arrangement. 

Direct calculations of these diagrams using standard rules of the diagram technique 
yield diverging expressions as a result of the integration over large momenta in subdiagrams 
containing Cooper loops (loops formed from two lines in the same direction). For example, 
we consider the first third-order diagram in Fig. |^ together with the corresponding diagram 
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Figure 3: Skeleton diagrams of the third order of perturbation theory for the irreducible 
vertex F 

in which is replaced with —p3. In expanded form they correspond to the sum of the 
diagrams in Fig. |^ where the dashed line is the interparticle interaction. It is easily to 
check that for short-ranged potential the first three diagrams cancel out leaving only the 
fourth diagram which contains a subdiagram diverging at large momenta, corresponding 
to a Cooper loop between two parallel dashed lines. However, this subdiagram is the first 
correction of a ladder series to the one of the dashed lines on the fourth diagram in Fig. ^ 
which gives a contribution to in the second order with respect to A. In the same way, the 
second diagram in Fig. ^ together with the corresponding diagram in which is replaced 
with — p3 in the sum is the correction to the second dashed line on the same diagram in 
Fig. |. 

These corrections only differ from the first term of the Born series for the scattering 
amplitude in that they contain the single-particle Green's functions in the medium G 
and not in vacuum G^'^\ However, at large momenta the difference between G and G^'^^ 
disappears so that the divergence in the diagram in Fig. |^ can be eliminated by changing 
from the bare interaction g to the scattering length a (renormalization procedure). This 
scattering length a is determined by the scattering amplitude of two particles in vacuum 
in the limit where the energies of the colliding particles tend to zero and can be obtained 
from the equation 




(13) 
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Figure 4: Skeleton "nonoriented" diagrams of the fourth order of perturbation theory for 
the irreducible vertex F. 
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After integrating in diagram Fig. ^ over the intermediate frequency of the Cooper 
loop UJ, we obtain the expression 



d^q f d^p 



1 - ^(6) - ^(6) 



(27r)4 J (27r)3 - (^i + 6) + ^5(sig<i + sign^) 

1 



{il- ^3 + i6sign^3){il - U + «5sign^4) 



(14) 



where = ^(p + ^), 6 = ^(-p + ^), ^3 = ^(q - s), ^4 = ^(q + s), pi = s + w, 
P3 = s — w. The integral over the internal momentum p of the expression in brackets 
diverges at large momentum. As we have already noted, this divergence is of the same 
type as the Born correction to the scattering amplitude. Therefore, the divergence in the 
considered third order diagram can be eliminated by replacing the bare interaction constant 
g with the scattering length a in the second order diagram in Fig. |^. 

To within the required accuracy the relationship between g and a can be obtained from 



(p!3|) and has the form 



47r Mtt 
— a + — a 

m \m 



Air 



m 



-a + 



47r ^2 
— a 

m . 



d^pduj 
d^p 



G(°)(cu,p)GW( 

1 



-p) 



(15) 



i2TiY 2e{p) + i5' 
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p2 

where e(p) = — . The renormahzation procedure is shown schematically as follows: 
2m 

Ma^^ + 2g^{GG)Il} ^ A'n + 2X^[{GG) - (G^^^G^"))]!!, (16) 



where (GG) corresponds to the first cofactor in formula (0), H corresponds to the second 
cofactor, 

= / pG^'\uj, -p)G(°)(^, -P) ^ 



and the factor 2 in second term in (16) results from the the contribution of the second 
diagram in Fig. |^. As we can easily see, the expression for 11, being integrated over 
frequency, gives the polarization operator 11 which appears in the first term in formula 
(16). The last term in formula (16) can be explicitly written as 

r(-) = _272^y/^/^ I .7) 

[mpFj J (27r)4i {27i)^n-^s + t5sign^^){n-^i + i5sign^4) 



+ 6) + «5(signei + sign6) 2e{p) + i6 ) ' 

This expression contains no divergences and can be integrated numerically. It can 
be checked that all other third-order diagrams contain no divergences and, as a result of 
numerical calculations, we obtain the final result for the third-order contribution to the 
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Figure 6: First third-order diagram from Fig. 3 showing diagram corresponding to the 
substitution p3 — > — Ps, in expanded representation. 



p-wave harmonic of the irreducible vertex: 

vpVf = -0.33AI (18) 

It should be noted that formula (0) contains no contribution from Hartree-Fock self-energy 
components in the second-order diagrams since this contribution corresponds to renormal- 
ization of the chemical potential. We also note that the appearance of a large numerical 
coefficient 0.33 (compared with the coefficient of 0.077 for the second-order contribution) 
is associated with the stronger angular dependence of the third-order diagrams (see Fig. 

This dependence is mainly determined by the first two diagrams in Fig. ^ and can be 
attributed to the existence of subdiagrams with Cooper loops. 

All divergences in the fourth-order diagrams can be eliminated in exactly the same 
way. For this purpose, in the third-order diagrams we need to express g in terms of a in 
accordance with formula (p!5D and in the second-order diagrams, express g in terms of a 
allowing for the term ~ a^, which can easily be obtained from equation (|13]). (This term 
is required to eliminate the divergences in the second diagram in Fig. ^.) As a result, the 
contribution of the fourth-order diagrams in Fig. ^ is given by 

upti"^^ = -0.39X\ (19) 

In order to calculate Fi to within A'* we also need to allow for the contribution of the 

~ (2) 

self-energy insertions of the second order in A in the second-order diagrams F^ , see Fig. 2. 
These contributions can no longer be reduced to renormalization of the chemical potential. 
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Figure 7: Diagram of the third order of perturbation theory containing a Cooper loop as 
sub diagram. 



They also result in the appearance of a non-trivial Z-factor and the effective mass m* [0: 

m* 2 

Z = l-XHn2- — = 1 + — (71n2- 1)A2 (20) 
m 15 

in the singular part of the single-particle Green's function which now also contains a regular 
part proportional to A^. By means of direct numerical calculations of the corresponding 
diagrams we can establish that the contribution of the latter to Fi is negligible. Thus, 
we finally obtain the following expression with the required accuracy in terms of A for the 
irreducible vertex in a Cooper channel with orbital momentum i = 1: 

m m \ m 

= -O.OTTA^ - 0.33A^ - 0.26A1 (211 



4 The retardation effects 

In order to determine the critical temperature in Section 2, in equation (1) we replaced 
the irreducible vertex Fi, which is a function of the incoming and outgoing frequencies and 
absolute values of momenta Ti{uji, pj), by its value at zero frequencies and momenta lying 
on the Fermi surface Ti{uji = 0, pj = pp). In this section we shall show that allowance for 
the difference between ri(ci;j,pj) and Ti{u!i = 0, pj = pp) (retardation effects) introduces 
a correction of the order of A'^ to the vertex Fi. In other words, these effects influence the 
numerical coefficient in the preexponential factor. 

Retardation effects are most conveniently taken into account using a method proposed 
in ip, . 2]. Omitting the appropriate procedures, which are a trivial generalization of the 
derivation of to the case of p-wave pairing, we arrive at the following integral equation: 

MO = -J de'^^^R,{^,am'), (22) 
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Figure 8: Dependence of the irreducible vertex in the second and third orders on the 
angle between the incoming and outgoing momenta ^1-3: T^'^^X'^ - solid curve, T^^^/X'^ - 
dashed curve. 



for which the condition for existence of a nontrivial solution determines the critical tem- 
perature Tel. The unknown function $1 in ([2^ ) can be related to the spectral density of 
the anomalous Green's function (or more accurately to its first harmonic in the expansion 
in terms of Legendre polynomials), and the kernel is given by 



m 



P ,f f'l "Y" 1 dEa(E,m) AQ + k^O'g' 

Ip(0-pK')I 

where (j{E, q) is related to r(ci; = uji — uj^, q = Pi — Ps) by 

r(z^„.,p) = y ^ , (24) 

" 

and the factor (p^ + fc^ — q^)/2pk is precisely the cosine of the angle between the incoming 
and outgoing momenta which picks up the first harmonic in the expansion (^) in terms of 
Legendre polynomials. 

Dividing the region of integration over ^' in equation (|2^) into three parts: |^'| < zep-, 
—Ep < ^' < —zefi and S,' > zep, where z is an arbitrary number satisfying the condition 
Tci <ti ZEp Ep, and integrating by parts (where the dependence on ^' in R{^,^') and 
$i(^') can be neglected in the first region and the hyperbolic tangent in the second and 
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third regions can be replaced by =f1, respectively), equation ( p2D can be reduced to the 
form 

m) = - In (^) $i(0)i?i(e,0) + lJdClr^(^)^ n<^>i(0) • (25) 



-ep 



(As was to be expected, the arbitrary constant z was dropped from this equation.) We 
introduce the new variable 

xiO = ^'^%T^ , (26) 



$i(0) In 



cl 



27£f 



which allows us to write the expression for the critical temperature in the form 



2e^ 

Tci = ep exp 



TT 



x(o). 



(27) 



where the function x(0 satisfies 

xio = Riim + \] ^e'ln ^ mi,Ox{0) 



(2J 



Since the kernel Ri contains the small parameter {R ~ A^), equation ( p2D can be solved 
by an iterative method. In the zeroth approximation we set: 



x^"He) = i?i(e,o). 

The first correction x*^^) is given by the integral on the right-hand side of (p8D with x = x'"'^^'- 



X 



(1) 



d^' In 



-ep 



Ef) di' 



;(i?l(e,O^l(e',0)) 



(29) 



and, as can easily be seen, begins with terms of the order A^. The leading term with respect 
to A in Xi is obtained if only the leading (~ A^) terms are retained in the kernel Ri in 
formula (|29|) . In this case, the spectral function a{E, q) is the same as the imaginary part 
of the operator 



a(E,q) 



-lmR{E, q) 
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Figure 9: Dependence of Td/ep on the gas parameter A. 

where q = q/pr, E = Em/p'jp. Direct calculations using formula (p9|) give x^^^ ~ 0.004A^ 
which is equivalent to adding AFi in formula (|^), 

ufATi ^ 0.004A^ (30) 

Note that a similar estimate of the contributions of retardation effects was made in 
where the authors used the step function approximation for the frequency dependence of 
the polarization operator. 



5 Critical temperature 



cl 

Collecting above results together [formulas (|21| ) and (^)], we obtain the following expres- 
sion for the critical temperature which is determined numerically to within two decimal 
places: 

Td = -e^ej.exp{-(0.077A2 + 0.33A2 + 0.26A^)"^} — > 
2 c / 13.0 42.0 I 

where the omitted terms have the order A. This last formula assumes A < 0.23 since for 
A = 0.23 the second- and third-order terms with respect to A in the exponential function 
in (|3T|) are the same. For 0.23 < A < 1 the fourth-order term with respect to A in ( ^T]) is 
smaller than the first two so that (|3T|) can accurately be rewritten in the form 

^ 2 r { 13.0 42.0 1 , , 

Tci ^ -e^sp exp <^ - — + — \ 32 

TT ^ ^1 A2 1 + 4.3A l + 4.3A2j ^ ' 
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This formula may be considered as an extrapolation of the expression for the critical 
temperature from A ^ 1 [formula (|3l|) ] to the region A < 1 [formula (|32D ]. The dependence 
Tci(A) is shown in Fig. 



6 Influence of bare p-wave scattering 

So far we have only considered s-wave scattering between particles, assuming that the 
interparticle potential is a point-like. However, as we have already pointed out, for the 
potential with a finite range, the problem will always contain scattering having an arbitrary 
orbital momentum £ whose amplitude for particles having momenta equal to the Fermi 
momentum pp may be estimated as /; ~ a{appY^ [^. From this it follows that with the 
required accuracy we can confine our analysis top- wave scattering {^ = 1). In this case, only 
two contributions will be important: a contribution of the order A^ fromp-wave scattering 
at the bare interparticle potential and a contribution of the order A'' corresponding to the 
diagram in Fig. |^d where one of the dashed lines corresponds to s-wave scattering and the 
other to p-wave scattering. More precisely, if the amplitude of p-wave scattering of two 
particles having momenta is written in the form 

/i = aia ^ — = aioA^, (33) 

ai is a numerical coefficient of order unity, the contribution of triplet scattering to the 
irreducible vertex Fi has the form 

up^i = aiA^(l + 0.008A). 

We can see that the fourth-order contribution with respect to A can be neglected and 
consequently the critical temperature is given by 

^ 2 ^ f 13.0 42.0 1 , , 

= '^'"P i-A^[l + (4.3 + aOA] + [1 + (4.3 + a,)AF / " ^''^ 
Nevertheless, we can specify a physical situation when the contribution of p-wave scat- 
tering can be neglected. This corresponds to the case when a shallow level having the 
orbital momentum £ = (resonance scattering) exists in the potential. In this case, the 
p-wave scattering amplitude is estimated as /i ~ to('"oPf)^5 where tq is the radius of action 
of the potential while the s-wave scattering length is given by a = (l/2m|£'|)^/^ ^ ro, 
where E is the discrete level energy (we assume that the condition \E\ ^ Sp is satisfied so 
that app -C 1). Then for ai in formula (^) we obtain the estimate 

and if the condition 

tti < A 

is satisfied, the contribution of the p-wave harmonic of the bare interparticle interaction 
can be neglected compared with the fourth order of the effective interaction which allows 
only for s-wave scattering. 
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7 Critical temperature in a magnetic field 



In this section we study the influence of an external magnetic field on the irreducible 
vertex Fi and consequently on the critical temperature Td to within terms of the order 
A^. As was shown in 



m 



in the leading approximation with respect to A in the model 
being studied the critical p-wave pairing temperature may increase appreciably if a static 
magnetic field is applied to the system. This is because for conventional singlet pairing 
the role of a magnetic field is always destructive due to the paramagnetic suppression of 
Cooper pairing caused by the flipping of one of the pair spins. However for triplet p-wave 
pairing no paramagnetic effect occurs so that the role of the magnetic field is not clear a 
priori. 

In our approach the mechanism for variation of T^i in a magnetic field is based on 
the magnetic field dependence of the many-particle effects which determine the effective 
interaction. On the one hand, as a result of a difference in the number of particles (and 
consequently Fermi momenta) with spins directed parallel and antiparallel to the field, 
the Kohn singularity increases sharply, causing an increase in Fi. On the other hand, 
i^Fi^i value decreases with increasing magnetic field because of a monotonic decrease in the 
number of particles with spin antiparallel to the field. (We recall that s-wave scattering can 
only occur between Fermi particles having different spin projections.) Competition between 
these two effects leads to an sharply non-monotonic dependence of the critical temperature 
on the magnetic field (more accurately, on the degree of polarization a = (n|— ra|)/(n|+n|)) 
with pronounced increase in T^i for small a, a maximum at intermediate a, and a decrease 
for a — >■ 1. (In this case of a completely polarized Fermi gas only p-wave scattering between 
parallel spins can only take place.) The dependence of F on the polarization a to within 

[0: 



second-order terms was calculated in 



f(2)(5) = -A' 



,21n2 -11 



^3 Vl + 



2/3 



1 + 



6 



3(2 In 2 



(35) 



where 



{6 + 1) 101n(5 + l) -52-3 



-{6^ + 26^ + 86 + 4) In 



6 



+ 



6-16 



In 



(5+1) 



Pfi \l-aj 

In the third order with respect to A the result can only be obtained numerically. The 
corresponding contribution is given by the diagrams in Fig. ^ where the spins on the 
outer lines are directed parallel to the field and those in the inner loops can be oriented 
either parallel or antiparallel to the field. The calculations (including renormalization of 
the diverging diagrams) are exactly the same as in the absence of a magnetic field and the 
result is shown in Fig. |10| (solid curve) which also contains the second-order contribution 
(P5|) for comparison (dashed curve). 
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a 

Figure 10: Dependence of the second- and third-order contributions to the irreducible 
vertex on the degree of polarization a: T^'^yX^ - dashed curve, F^^^A^ - solid curve. 

It can be seen that the maximum of T\ (a) is obtained at = 0.48 whereas T\ (a) 
decreases monotonically. Thus, the maximum of is determined by competition between 
increasing F^ (a) and decreasing F^^ (a). For typical A this is in the region of a ~ 0.4. 
Graphs of the critical temperature as a function of the degree of polarization are shown 
in Fig. |ll] for typical values of A. For A = 0.6 the value of Td at the maximum is 
approximately six times large than the value of Td in the absence of the field. In the latter 
case the maximum is mainly determined by the second order and is reached at A ~ 0.45. 

8 Discussion of the results 

The experimental search for nontrivial pairing with / 7^ gl in isotropic Fermi systems has 
recently been actively pursued. Until recently the main candidate was ^He-^He mixture. 
So far superfluidity has not yet been observed in this system although temperatures of 
the order of 97 fiK have been achieved experimentally ||Tl|. In the concentration range 
X < Xo ~ 3% the scattering length in the mixture correspond to the attraction so that 
singlet s-wave pairing may be achieved. The critical temperature is given by formula (^) 
allowing for 

Ef = SfqX^'^, Pf = PfoX^^^, 

where Efo, Pfo are the Fermi energy and Fermi momentum of pure ^He. According to 
estimates made in we have 

maxTeo = T,o(l%) ^ lO'^K. 

The authors of [13] predict an even lower critical temperature: 

maxT,o = T,o(2%) ^ 4 ■ lO"*' - lO'^K. 
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Figure 11: Dependence of Tci/sf on the degree of polarization a for various A. 



13 



as 



Note that the value T^o ~ 10~^ was extracted from spin-diffusion experiments in 
a fitting parameter to describe magnetostriction experiments and T^q ~ 4 ■ 10~^K was 
obtained in spin diffusion experiments. It should be noted that for a given concentration x 
the gas parameter of the theory appox^^^ depends weakly on pressure. Hence the pressure 
cannot be considered as an instrument to obtain optimum parameters for s-wave pairing. 

For high concentrations (x > Xq) the scattering length changes sign a > and s- 
wave pairing becomes impossible. Nevertheless, in this case the subsystem of ^He atoms 
may become superfluid, but now with respect to p-wave pairing. The critical temperature 
is given by formula (0) with A replaced by Ax^/^ and ep replaced by epox^^'^- It has a 
maximum at P = 10 atm when the maximum ^He concentration of 9.5% is achieved. Figure 

gives the dependence of T^i on the concentration calculated by using the extrapolation 
formula (0). At maximum concentration x = 9.5% the temperature is of the order of 
10~^. A further increase in in solution may occur in strong magnetic field. For example, 
at X = 9.5% the maximum of Td in a field is more than six times that in the absence of a 
field, leading us to experimentally measurable temperatures of 6 • 10^^. 

Recently the properties of trapped Bose-condensed gases of alkali elements (^^Na, 'Li, 
^^Rb) have been studied intensively. A combination of laser and evaporative cooling in 
magnetic traps can reach gas-phase densities of the order of 10^^ — 10^^ cm~^ and temper- 
atures of the order of 10^^ — 10^^ . In addition these elements may have a anomalously 
large scattering length a of quasi-resonant origin. For Rb and Na the scattering lengths 
are positive. It is also found that the scattering length may cover a broad spectrum of 
values from negative to positive as a result of the Feshbach effect. This effect was observed 
for 23Na [0. 



A logical continuation of studies of Bose condensation in trapped gases of alkali atoms 
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Figure 12: Dependence of Td on the concentration x in a ^He-^He mixture for various 
degrees of polarization: a = 0.2 - solid curve, a = 0.1 - dashed curve, and a = 0.0 - 
dotted curve. 

would be to obtain superfluidity in low density fermionic system in restricted geometry. 
The case of a negative scattering length makes it possible to achieve s-wave pairing with a 
transition temperature determined by formula (H). For ^Li, for example, we have a = —2.3- 
10~^A < 0. Thus, for n ~ lO^^cm"^ the critical temperature Tco is of the order of ~ 10~^K. 
Note that because of the Pauli principle the wave function of an s-wave Cooper pair should 
be antisymmetric with respect to the interchange of quantum numbers characterizing the 
internal state of the atoms forming the pair. These numbers are indices determining the 
multiplet component of the hyperfine interaction for the case of zero field (optical trap) or 
weak magnetic field. They can also correspond to the projections of the nuclear spin when 
the strong external magnetic field of the trap destroys the hyperfine coupling. Thus, s-wave 
pairing can only take place between atoms of different gas components. This imposes a 
very stringent constraint on the closeness of their densities. From the experimental point 
of view we should have: \ni — n2|/(^i + ^2) < Tco/sp <^ 1 In the opposite case the Cooper 
pair would have a velocity higher than the critical velocity Vc ~ Tcq/pf-, and hence will be 
destroyed. From the experimental point of view, it may therefore prove difficult to achieve 
this type of pairing experimentally. For p-wave pairing, a Cooper pair may be formed by 
atoms of the same component (an analog of the A2 phase in superfiuid ^He). Note that 
the superfiuid transition temperature in the triplet case may be increased substantiafiy 
by utilizing the presence of several components in the trap. This increase is similar to 
the increase in Td in a magnetic field and is associated with the idea of the separation of 
channels: Cooper pairing is achieved between particles of one component as a result of the 
polarization of the other components. In this case, it is possible to obtain a superfiuid p- 
wave pairing temperature of the order of 10^'' — 10^^, which is quite feasible experimentally. 
By virtue of this fact this type of pairing may be quite promising from the experimental 
point of view. 
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